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Introduction 


Recent  low-altitude  satellite  measurements  of  electrons  in  the  drift  loss  cone 

between  L  ~  1.6  and  L  ~  1.8  display  several  characteristics  which  suggest  that  the 

electrons  were  precipitated  by  a  resonant  interaction  with  waves  from  a  ground-based 

VLF  transmitter  [Vampola  and  Kuck,  1978;  Imhof  et  al.  1974a,  b].  The  easternmost  edge 

of  the  pitch-angle  scattering  region  can  be  determined  from  the  size  of  the  observed  loss 

cone  in  the  particle  pitch-angle  distribution  under  the  following  assumptions  (1)  the 

pitch-angle  scattering  is  sufficiently  strong  to  fill  the  drift  loss  cone  to  at  least  the  edge 

of  the  local  bounce  loss  cone,  (2)  the  atmosphere  removes  all  particles  which  have  a 

mirror  altitude  below  100  km  at  the  location  where  they  have  been  pitch-angle  scattered 

down  from  higher  mirror  altitudes  and  (3)  that  no  further  scattering  occurs  as  the 

particles  drift  eastward.  Using  the  observed  value  of  the  loss-cone  angle  of  electrons  in 

the  drift  loss  cone,  Vampola  and  Kuck  [1978]  showed  that  the  particles  could  be  mapped 

westward  in  longitude  to  the  region  in  which  they  were  precipitated  (see  Vampola  and 

Kuck  [1978]  for  a  discussion  of  this  mapping  procedure  and  the  limits  of  accuracy).  For 

precipitation  events  in  the  L  range  from  1.6  to  1.85,  they  found  that  most  of  the  events 

mapped  back  to  a  small  region  of  longitude  near  60°  E.  Assuming  a  density  model  with  a 
-4 

L  dependence  and  ducted  propagation,  Imhof  et  al.  r1974a,  b]  found  that  the  particles 
they  observed  could  resonate  with  ducted  waves  at  a  frequency  of  10  kHz,  which  is  very 
close  to  the  lowest  frequency  used  by  the  world-wide  Omega  navigation  transmitters  at 
10.2  kHz.  However  the  authors  did  not  make  a  positive  identification  of  the  source  of 
the  waves,  nor  did  they  show  that  such  ducted  waves  occurred  in  this  region  of  the 
plasmasphere. 

The  measurements  reported  by  Vampola  and  Kuck  [1978]  and  by  Imhof  et  al. 
[1974a,  b]  have  the  following  characteristics  in  common:  all  of  the  observations  exhibit 
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a  very  narrow  energy  distribution  and  the  energy  decreases  with  increasing  L-shell.  An 
example  from  the  OV1-19  data  is  shown  in  Fig.  1.  Furthermore,  Vampola  and  Kuck 
[1978]  using  high-resolution  pitch-angle  data  show  that  the  particles  entered  the  drift 
loss  cone  in  a  narrow  longitude  range  which  was  typically  just  to  the  west  of  60°  E.  The 
first  and  second  characteristics  suggest  that  there  was  a  resonant  interaction  with  a 
monochromatic  signal.  The  third  characteristic  indicates  that  the  signal  was  generally 
confined  to  one  locality  on  the  ground,  presumably  a  VLF  transmitter  in  the  Soviet 
Union.  In  this  paper  we  investigate  the  resonant  interaction  with  a  ray-tracing  analysis 
and  use  data  from  a  ground-based  VLF  receiver  to  show  that  the  sources  of  the  waves 
were  the  VLF  transmitter  UMS  located  near  Gorki,  USSR  and  the  NWC  transmitter 
located  at  North  West  Cape,  Australia. 

Figure  2  shows  a  diagram  of  the  interaction  mechanism  in  the  inner  magneto¬ 
sphere.  A  ground-based  VLF  transmitter  at  a  midlatitude  location  excites  VLF  waves 
which  propagate  worldwide  in  the  earth-ionosphere  waveguide.  Some  of  the  VLF  energy 
couples  into  the  ionosphere  and  enters  the  magnetosphere  where  both  ducted  and 
nonducted  propagation  takes  place.  The  maximum  intensities  in  the  magnetosphere 
occur  near  the  longitude  of  the  VLF  transmitter.  The  nonducted  rays  entering  between 
the  magnetic  dipole  latitudes  of  23  and  31°  cross  the  magnetic  equator  between  L  ~ 
1.5  and  L  ~  1.8  with  wave-normal  angles  near  63°  with  respect  to  the  magnetic-field 
direction.  In  a  region  about  the  equator  these  VLF  waves  can  interact  with  100  to  500 
keV  electrons  by  a  doppler-shifted  cyclotron  resonance.  Pitch-single  scattering  resulting 
from  the  resonance  process  will  cause  the  electrons  to  be  precipitated  into  either  the 
drift  loss  cone  or  the  local  bounce  loss  cone  where  they  can  be  observed  by  a  low- 
altitude  satellite  such  as  OV1-19. 
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Fig.  1.  Electron  Flux  vs.  Time  plots  of  data  obtained  in  the  drift  loss  cone  by  the 
OV1-19  satellite  in  1969.  The  location  of  the  peak  in  L  increases  monotoni- 
cally  with  decreasing  electron  energy.  The  scatter  in  data  points  is  due  to 
pitch-angle  sampling. 
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Fig.  2. 


Cartoon  presenting  the  interaction  of  coherent  waves  from  a  ground  based 
transmitter  with  electrons  near  the  geomagnetic  equator.  Waves  from  the 
transmitter  propagate  southward  in  the  ionospheric  cavity  and  leak  up 
through  the  ionosphere  to  the  equatorial  region  where  they  interact  with 
northward  travelling  electrons.  The  electrons  lose  perpendicular  momentum 
to  the  waves,  thereby  lowering  their  pitch-angles.  Some  are  lost  into  the 
atmosphere;  the  remaining  electrons  then  have  a  pitch-angle  distribution 
which  is  characteristic  of  the  location  at  which  the  atmospheric  interaction 
occurred.  The  distribution  is  then  observed  by  a  low  altitude  satellite,  such 
as  OV1-19,  as  the  electrons  drift  eastward  to  the  South  Atlantic  Anomaly 
(where  the  distribution  is  changed  back  to  that  characteristic  of  the  SAA). 
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Ray  Tracing  Analysis 

There  are  several  reasons  why  nonducted  propagation  is  assumed  instead  of  ducted 
propagation.  Examining  the  doppler  shifts  of  waves  from  a  VLF  transmitter  arriving  in 
the  topside  ionosphere  of  the  conjugate  hemisphere  Cerisier  [1974]  and  Edgar  [1976] 
found  that  below  L  -  2  all  satellite  observations  can  be  satisfactorily  explained  by 
nonducted  propagation.  Ground-based  observations  of  low  L-shell  whistlers  are  highly 
dependent  upon  magnetic  storms  for  the  production  of  density  gradients  necessary  for 
the  proper  ducting  of  the  ray  paths  [Singh,  1976].  All  of  the  observations  reported  by 
Vampola  and  Kuck  [1978]  occurred  during  magnetically  quiet  or  moderately  disturbed  (K^ 

<  3)  conditions  under  which  the  density  structures  producing  whistler-mode  ducts  are 
not  likely  to  exist.  Thus,  nonducted  VLF  propagation  is  the  assumed  mode  for  L  <  2. 

Although  observed  far  to  the  east  of  the  transmitter,  the  electrons  were 
presumably  precipitated  near  the  longitude  of  the  transmitter  wnere  its  signal  strength 
in  the  magnetosphere  would  be  strongest.  The  electrons  then  drift  eastward  at  a  rate 
which  depends  upon  their  energy,  equatorial  pitch  angle  and  L-shell.  The  electrons 
observed  by  the  satellite  in  a  very  short  span  of  time  (typically  ~  30  sec  for  the  data  as 
shown  in  Fig.  1)  were  in  fact  precipitated  more  than  an  hour  earlier  with  the  lower- 
energy  particles  precipitating  tens  of  minutes  before  the  higher-energy  particles. 
Detailed  calculations  are  presented  below. 

VLF  ray-tracing  calculations  were  performed  for  a  model  plasmasphere  in 
diffusive  equilibrium  [Angerami  and  Thomas,  196 4].  The  diffusive-equilibrium  model  used 
a  90  to  10  percent  H+  to  0+  mixture  at  1600°  K  at  1000  km.  The  model  is  parameterized 
by  the  electron  density  at  1000  km.  The  ray-tracing  calculations  were  done  for  densities 
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at  1000  km  of  5  x  10^,  104,  and  2  x  104  em*^.  The  solid  lines  in  Fig.  3  show  the  L-shell 
profile  of  the  electron  density  at  the  magnetic  equator  for  these  models. 

Assuming  that  the  interaction  between  the  particles  and  the  nonducted  waves 
occurs  at  the  equator,  the  important  parameter  for  the  resonance  calculation  determined 
from  the  ray-tracing  analysis  is  the  wave-normal  angle  of  the  whistler-mode  wave  in  the 
interacting  region.  The  wave-normal  angle  at  the  equator  determined  for  the  three 
models  is  shown  in  Fig.  4.  The  wave-normal  angle  lies  between  62  and  64°  and  is 
essentially  independent  of  the  model  density  over  most  of  the  L-shell  range  of  interest. 

Resonance  Calculations 

An  electromagnetic  wave  propagating  in  the  whistler  mode  resonates  with 
electrons  in  the  inner  zone  through  the  doppler-shifted  cyclotron  resonance.  The  general 
resonance  condition  is 


w  -  k  ■  v  =  mn  (1) 

where  u>  is  wave  frequency,  k  is  the  wave  vector,  v  is  the  electron  velocity,  m  is  the 
order  of  the  resonance,  and  ft  is  the  electron  cyclotron  frequency.  For  the  first-order 
cyclotron  resonance  for  relativistic  electrons  this  equation  may  be  rewritten 

w  (1  +  0  ji  cos  e  cos  a)  =  ft  (1-0  (2) 

where  0  =  v/c,  a  is  the  electron  pitch  angle,  0  is  the  wave  normal  angle,  and  m  is  the 
index  of  refraction  for  the  whistler  mode. 


10 


Fig.  3. 


Equatorial  electron  density  as  a  function  of  L  calculated  using  a  diffusive- 
equilibrium  model  and  the  base  densities  shown.  Electron  densities  for  two 
electron  precipitation  events  are  plotted  as  squares  and  circles.  See  text  for 
details. 
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4.  Wave-normal  angle  at  the  equator,  assuming  unducted  propagation,  for  the 
three  models  of  Fig.  3.  The  angle  is  essentially  independent  of  the  base 
density  assumed. 
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The  index  of  refraction,  which  is  the  ratio  of  the  velocity  of  light  to  the  phase 
velocity  of  the  wave,  may  be  written 

J  1/2 

4ff0-c&"e  rgl  (3) 

2  2 

where  cjp  =  Ne  /m0E0-  The  electron  plasma  frequency  depends  only  on  the  electron 
density,  N,  in  the  diffusive  equilibrium  model  of  the  plasmasphere. 

For  an  interaction  at  the  magnetic  equator  with  a  specific  frequency  f  =  gj/2  ir  on 
a  given  L-shell  which  determines  C2q,  with  waves  of  wave-normal  angle  e  which  is 
determined  from  the  ray-tracing  calculations  and  is  a  function  of  L,  Eq.  2  is  a 
relationship  between  the  particle  energy  which  is  determined  by  (J ,  its  pitch  angle  a  and 
the  electron  density  through 

The  equatorial  pitch  angle  used  in  the  resonance  calculation  is  the  pitch  angle  at 
the  edge  of  the  stable-trapping  region.  This  represents  particles  which  have  just 
survived  an  azimuthal  drift  through  the  South  Atlantic  Anomaly.  As  shown  in  Fig.  5  the 
pitch  angle  that  defines  the  limit  to  the  stable-trapping  region  is  a  function  of  L. 

We  are  finally  left  with  a  relationship  between  the  equatorial  electron  density  as 
a  function  of  L  and  the  energy  of  the  corresponding  resonant  electrons. 

Since  the  electrons  are  observed  in  our  experiment  we  use  Eq.  2  to  compute  the 
electron  density  as  a  function  of  L.  The  L-shell  at  which  the  particle  count  rate  reaches 
a  peak  in  each  energy  channel  is  taken  to  be  the  L-shell  of  the  interaction  for  that 
energy.  The  L-shell  dependence  of  the  electron  energy  extracted  from  the  data  shown  in 
Fig.  1  is  shown  in  Fig.  6. 
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Fig.  5.  Equatorial  pitch-angle  as  a  function  of  L  for  particles  in  the  local  bounce  loss 

cone  and  stable  trapping  zone.  Particles  with  intermediate  equatorial  pitch- 
angles,  unhatched  area,  will  mirror  above  the  atmosphere  over  some  portion 

of  longitude  but  will  encounter  the  atmosphere  prior  to  completing  360°  of 


drift. 


The  electron  density  calculated  from  Eq.  2  using  the  procedure  described  above  is 

plotted  as  circles  on  the  graph  in  Fig.  3.  The  data  generally  follow  the  curve  which  is 

4  -3 

the  model  for  diffusive  equilibrium  with  a  base-level  electron  density  of  10  cm  at 
1000  km. 

The  L-shell  dependence  of  the  peak  count  rate  for  each  energy  channel  for  all  of 

the  events  detected  by  satellite  OVl-19  is  shown  in  Fig.  7.  The  predicted  energy 

dependence  for  diffusive-equilibrium  models  with  base-level  electron  densities  of  104 
4  -3 

and  2  x  10  cm  at  1000  km  altitude  are  drawn  as  dashed  lines  in  Fig.  7. 

The  slope  of  the  data  is  in  excellent  agreement  with  the  slope  of  the  predicted 
dependence  and  most  of  the  events  lie  within  the  limits  defined  by  these  density  limits. 
Departures  from  diffusive  equilibrium  lead  to  changes  in  the  slope  of  the  energy 
dependence  vs  L-shell  curve. 

Well  defined  peaks  in  the  data  can  be  used  to  estimate  the  electron  density  profile 
by  the  technique  used  above  to  obtain  the  data  points  shown  on  Fig.  3. 

As  mentioned  earlier  Imhof  et  al.  [1974a,  b]  observed  electrons  in  the  drift  loss 
cone  with  an  energy  vs  L-shell  dependence  similar  to  those  observed  by  OVl-19.  We 
have  replotted  their  data  in  Fig.  8  with  bands  showing  the  predicted  energy  dependence 
for  diffusive  equilibrium  models  with  base  level  electron  densities  of  104  and  2  x  104 

_3 

cm  at  1000  km  for  nonducted  waves  between  15  kHz  and  18  kHz.  This  model  gives  a 
reasonable  fit  to  both  the  magnitude  and  slope  of  the  data. 

Since  their  data  was  taken  just  west  of  the  South  Atlantic  Anomaly,  the  particles 
could  have  been  precipitated  by  any  of  the  transmitters  that  are  active  in  the  VLF 
communications  band. 
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Data  similar  to  Fig.  6,  but  includes  all  of  the  events  identified  in  the  OVl-19 
data  set.  The  dashed  lines  are  the  result  of  calculations  using  a  diffusive- 
equilibrium  model  and  the  procedure  outlined  in  the  text 


1.45  1.5  1.55  1.6  1.65  1.7  1.75  1.8  1.85 

L 

Plot  similar  to  Fig.  7,  but  using  the  data  of  Imhof  et  a J.  (1974a,  b).  Solid  lines 
are  the  calculations  of  Imhof  et  al.  The  diffusive-equilibrium  model  with  the 
procedure  in  the  present  work  (grey  bands)  gives  a  better  fit  to  the  slope. 


Comparisons  of  VLF  Transmissions  and  Particle  Data 


Electrons  precipitated  into  the  drift  loss  cone  will  drift  eastward  until  they 
interact  with  the  atmosphere  and  are  lost,  usually  in  the  region  of  the  South  Atlantic 
Anomaly.  The  rate  at  which  they  drift  eastward  is  a  function  of  energy,  L-shell,  and 
equatorial  pitch  angle. 

The  drift  period,  2ir/^g,  can  be  calculated  using  the  formula  given  by  Schulz  and 
Lanzerotti  [1974]: 

2:r/n3  =  — (3L/2  iry)  ( y 2  -  1)  (c/a)2(m0/qBQ)  D(y)/T(y)  (4) 

where  y  is  the  relativistic  mass  factor,  c  is  the  velocity  of  light,  a  is  the  radius  of  the 
earth,  mQ  is  the  electron  rest  mass,  q  is  the  electron  charge,  BQ  is  the  equatorial 
magnetic  field  at  the  surface  of  the  earth,  and  y  is  sin  a.  D(y)  and  T(y)  are  given  by: 

T(y)  =  T(0)  -  1/2  [T(0)  -  T(l)]  (y  +  y1/2) 

D(y)  *  1 4T(0)  -  [3T(0)  -  5T(1)]  y  -  [T(0)  -  T(l)]  (y  In  y  +  y1/2)  |  /12 

T(0)  -  1.3802 
T(l)  a  0.7405 

The  effect  of  the  difference  in  drift  velocity  as  a  function  of  energy  and  L-shell 
will  be  to  produce  a  dispersion  in  time  for  the  various  energy  components  simultaneously 
precipitated  at  the  same  longitude.  Since  the  difference  in  drift  velocity  can  be  quite 
significant,  the  particles  observed  at  a  given  longitude  may  have  been  precipitated  at 


widely  varying  times.  As  an  example,  particles  observed  at  120°  E.  can  include  380  keV 
electrons  precipitated  22  minutes  earlier  (assuming  a  precipitation  at  44°  E.)  and  139 
keV  electrons  precipitated  50  minutes  previously.  In  spite  of  the  dispersion  which 
occurs,  it  is  possible  to  attempt  to  correlate  the  precipitation  with  the  operating 
schedule  of  the  assumed  transmitter,  UMS. 

Using  the  drift  period  obtained  from  Eq.  4  the  particles’  longitude  can  be  traced 
backward  in  time  from  the  longitude  of  the  observation.  The  drift  period  is  a  function  of 
the  particle  energy  through  y,  the  equatorial  pitch  angle  a  through  y,  and  the  L-shell. 
The  equatorial  pitch  angle  is  taken  to  be  the  equatorial  pitch  angle  of  electrons  which 
are  locally  mirroring  at  the  satellite  in  the  drift  loss  cone  at  the  location  of  the 
observation. 

VLF  recordings  from  Byrd  Station,  Antarctica  (private  communication,  R. 
Helliwell)  were  examined  for  time  periods  when  data  in  the  drift  loss  cone  were  available 
from  the  OV1-19  satellite.  There  were  three  periods  when  both  particle  data  and  VLF 
recordings  at  the  appropriate  times  were  available.  The  time  required  for  these 
electrons  to  drift  from  the  assumed  longitude  of  origin,  44°  E.,  to  the  location  of  the 
satellite  data,  140  to  180°  E.,  is  typically  30  minutes  to  an  hour,  the  exact  time 
depending  on  the  total  longitude  through  which  the  particles  drift.  Hence,  the  recording 
must  actually  cover  some  time  preceding  the  particle  observation  times.  The  three 
periods  occurred  on  August  14,  1969  at  about  19:30  UT  and  on  August  15,  1969  at  about 
19:45  and  22:20  UT. 

The  data  shown  in  Fig.  9  were  obtained  at  177°  E.  on  August  15,  1969  at  about 
19:45  UT.  The  L-shell  dependence  of  the  peaks  in  each  energy  channel  is  shown  in  Fig. 
10  and  the  equatorial  electron  density  obtained  from  the  resonance  analysis  is  plotted  as 
squares  in  Fig.  3. 
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Fig.  9. 
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Data  similar  to  Fig.  1,  for  August  15,  1969. 
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Fig.  10.  Plot  similar  to  Fig.  6,  but  for  the  data  of  Fig.  9. 
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The  drift  rate  corresponding  with  each  energy  channel  was  computed  from  Eq.  4 
for  the  peaks  as  identified  in  Fig.  10.  The  resulting  longitude  trace  as  a  function  of 
Universal  Time  is  plotted  in  Fig.  11  for  each  energy  channel.  The  traces  converge  on  the 
longitude  of  the  satellite,  177°  E.,  at  the  time  of  the  observation,  19:45  UT.  They  show 
the  longitude  of  the  detected  particles  as  a  function  of  time  prior  to  the  observation. 

A  VLF  spectrogram  covering  the  frequency  range  from  0  to  20  kHz  was  obtained 
from  tape-recorded  data  taken  at  the  Byrd  Antarctica  VLF  receiver  site  on  August  15, 
1969  during  the  time  period  shown  in  Fig.  11.  The  stations  identified  from  this  record 
are  shown  in  Table  1.  UMS  and  NDT  were  identified  by  Morse  Code  call  signs  in  the 
transmissions.  The  others  were  identified  from  their  known  transmission  frequencies 
during  August  1969.  The  only  stations  operating  below  20  kHz  at  that  time  and  located 
at  longitudes  between  the  South  Atlantic  Anomaly  east  to  the  satellite  location  were 
GBR,  JXZ,  UMS,  and  RPS.  The  signal  strengths  from  JXZ  and  GBR  at  Byrd  are  very 
weak  and  the  stations  are  located  well  to  the  west  of  the  longitude  to  which  the  particles 
are  traced  back. 

Thus  UMS  and  RPS  are  the  likely  candidates  to  account  for  the  observed  electron 
precipitation.  UMS  is  located  near  Gorki,  east  of  Moscow,  at  approximately  56.2°  N.  and 
44°  E.  [Watt,  1967%  The  location  of  RPS  is  unknown,  however  Aksenov  et  al.  [1970, 
1975]  report  using  a  transmitter  at  16.1  kHz  at  43°  51*  N.  and  77°  11'  E.  for  satellite 
experiments  conducted  in  1967  and  1969.  We  shall  assume  that  RPS  is  located  near  77° 
E.  and  see  if  this  is  consistent  with  the  OV1-19  particle  measurements. 

We  assume  that  the  longitude  width  over  which  the  interaction  occurs  is  ~  15° 
about  the  longitude  of  the  station.  Although  somewhat  arbitrary  it  is  based  upon  our 
analysis  of  ground-based  transmitter  signals  detected  by  a  VLF  receiver  aboard  the  S3-3 
satellite. 
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Fig.  11.  Traceback  in  time  and  longitude  of  electrons  observed  by  the  OVl-19 
satellite  at  about  19:45  UT  on  August  15,  1969,  at  177°  east  longitude.  The 
dots  represent  the  longitude  at  which  the  bounce  loss-cone  angle  matches  the 
observed  cutoff  in  pitch-angle  in  the  electron  distribution.  The  boxes  marked 
"UMS"  denote  the  times  during  which  UMS  was  transmitting.  They  are 
centered  at  the  presumed  longitude  cf  UMS  and  have  a  width  of  ±  15°  (the 
assumed  maximum  longitudinal  distance  over  which  UMS  is  effective  in 
precipitating  electrons). 
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Table  1 


VLF  Transmitters  Identified  in  Byrd,  Antarctica 
VLF  Data  on  August  15,  1969 


Frequency,  kHz 

Station 

16.0 

GBR 

16.2 

RPS  (UMS) 

16.4 

JXZ 

17.1 

UMS 

17.4 

NDT 

17.8 

NAA 

21.4 

NSS 

22.3 

NWC 

23.4 

NPM 

24.0 

NBA 

E.  Long. 

N.  Lat. 

356.5 

55.0 

13.0 

66.4 

44.0 

55.8 

139.5 

35.7 

292.7 

44.7 

283.6 

39.0 

114.2 

-22.8 

202.2 

21.4 

280.4 

9.1 
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The  operating  time  periods  obtained  from  the  VLF  spectrograms  are  plotted  on 
Fig.  11  as  boxes  whose  abscissae  cover  the  time  periods  of  the  transmissions  and  whose 
ordinates  are  ±15°  about  the  longitude  of  the  transmitters.  Between  18:32  and  18:45 
UT  the  transmissions  from  UMS  and  RPS  are  coincident  in  time. 

The  particles  detected  by  OV1-19  at  19:45  UT  must  have  interacted  with  the 
waves  from  the  transmitters  within  the  longitude-time  boxes  shown  in  Fig.  11.  The 
traceback  longitude  for  the  particles  in  this  model  would  be  the  longitude  at  which 
particles  in  a  given  energy  channel  exited  from  the  last  box  they  traversed.  Thus  from 
Fig.  11  the  traceback  longitude  for  the  139  keV  particles  would  be  87°  E.  while  all  of  the 
other  energies  would  trace  back  to  59°  E. 

We  have  been  concentrating  on  particle  data  obtained  while  the  OV1-19  was  in  the 
drtft  loss  cone  since  such  data  gives  a  clear  upper  limit  to  the  length  of  time  between 
the  wave-particle  interaction  responsible  for  lowering  the  electron  mirror  points  and  the 
time  of  observation.  However,  the  effects  of  these  wave-particle  interactions  are  also 
observable  in  data  obtained  at  higher  altitudes.  Figure  12  presents  data  obtained  at 
19:24  UT  just  prior  to  that  shown  in  Figs.  9  to  11.  Here  we  see  that  electrons  have  been 
pitch-angle  scattered  from  the  stable-trapping  region  to  lower  pitch  angles.  Pitch  angles 
from  the  stable-trapping  boundary  to  the  angle  at  which  electrons  will  mirror  at  100  km 
at  the  longitude  of  UMS  are  smoothly  filled  in.  At  lower  angles,  there  is  still  some  flux 
present  (including  fluxes  in  the  local  bounce  loss  cone)  but  the  slope  of  the  flux  vs.  pitch- 
angle  curve  abruptly  changes.  If  the  particles  were  scattered  in  the  vicinity  of  UMS.  this 
is  a  feature  we  would  expect.  The  drift  time  for  these  electrons  from  59°  E.  to  the 
longitude  of  observation  is  76  minutes.  UMS  was  transmitting  at  the  required  time, 
18:08  UT,  as  shown  in  Fig.  11. 
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Fig.  12.  Data  obtained  at  higher  altitude  just  prior  to  that  shown  in  Figs.  9-11.  Local 
pitch-angles  have  been  transformed  to  equatorial  pitch-angles  for  comparison 
with  the  local  bounce  loss  cone  (BLC),  the  bounce  loss  cone  at  the  longitude 
of  UMS  (UMS)  and  the  bounce  loss  cone  at  the  field  minimum  at  the  South 
Atlantic  Anomaly  (SAA).  Note  that  the  distribution  is  smooth  from  the 
stable  trapping  angle  to  the  "UMS"  angle  and  drops  abruptly  there,  as  would 
be  expected  if  particles  were  strongly  pitch-angle  scattered  at  the  longitude 
of  UMS  and  drifted  to  the  longitude  of  observation  (190°)  with  only  small 
perturbations  of  their  pitch-angles. 
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The  data  from  the  observation  on  August  14,  1969  is  shown  in  Fig.  13  in  the  form 
of  a  traceback  in  time  and  longitude  from  the  point  of  observation  of  the  particles  by  the 
OV1-19  to  the  time  and  longitude  of  their  last  interaction  with  the  atmosphere  (and 
presumed  longitude  of  origin).  On  August  14,  particles  were  observed  in  only  three 
energy  channels.  The  transmission  times  of  UMS,  as  determined  from  the  Byrd 
recordings,  are  shown  as  rectangular  boxes  as  in  Fig.  11.  All  three  channels  trace  back 
to  periods  of  transmissions.  The  250  keV  channel  had  a  poorly  defined  local  loss-cone 
with  a  shape  that  indicated  that  a  spread  in  longitude  was  contributing  particles  to  it. 
The  limits  of  this  spread,  as  determined  from  the  pitch-angle  distribution  and  the 
altitude  of  the  satellite  (particles  originating  west  of  48°  E.  would  have  been  mirroring 
above  the  satellite),  were  48°  and  58°  E.  The  48°  E.  location,  using  a  traceback  energy 
of  250  keV,  requires  that  the  channel  be  responding  to  electrons  of  energies  of  250  keV 
and  below.  Particles  with  higher  energy  would  have  already  drifted  past  the  location  of 
the  satellite.  The  250  keV  particles  would  have  been  precipitated  by  the  trailing  end  of 
the  third  transmission.  Particles  with  lower  energies  could  trace  back  to  either  the  same 
longitude  at  earlier  times  or  the  same  time  at  more  eastern  longitudes.  Particles  at  the 
nominal  lower  energy  limit  of  the  channel  response,  225  keV,  trace  back  to  the  58°  E. 
limit  at  the  same  time  as  the  250  keV  electrons  trace  back  to  48°  E.  Hence,  the  channel 
is  responding  to  particles  with  energies  in  the  lower  half  of  the  channel  response,  all  of 
which  came  from  the  trailing  end  of  the  transmission  and  a  spread  in  longitudes 
corresponding  to  the  variation  in  drift  speed.  The  poor  loss-cone  definition  is  the  result 
of  the  easternmost  edge  of  the  precipitation  region  not  being  represented  in  the  entire 
energy  spread  accepted  by  the  channel.  A  somewhat  similar  effect  occurs  in  the  192  keV 
channel.  Here  the  channel  is  responding  to  the  leading  edge  of  the  third  transmission  and 
the  response  is  primarily  to  higher  energy  electrons  which  have  drifted  faster  than  the 
nominal.  The  lower  energy  particles  have  not  yet  arrived  at  the  location  of  the  satellite. 
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UNIVERSAL  TIME,  hrmin 

Fig.  13.  Presentation  similar  to  Fig.  11,  but  for  data  obtained  on  August  14,  1969. 
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The  channel  is  then  slightly  biased  to  longitudes  to  the  west  of  the  easternmost  edge  of 
the  interaction  region. 

The  same  effect,  i.e.,  a  channel  responding  to  a  particle  distribution  which  is  not 
uniformly  represented  across  the  entire  energy  response  of  the  channel,  could  be  the 
cause  of  the  anomalous  points  on  Figs.  3  and  10.  If  the  139  keV  channel  is  responding 
primarily  to  the  leading  edge  of  the  second  transmission  sequence  (starting  at  18:32  in 
Fig.  11),  an  effective  energy  of  160  keV  results  from  Fig.  11.  This  is  within  the  energy 
band  of  that  channel.  The  resulting  calculated  points  for  Figs.  3  and  10  then  become 
consistent  with  the  rest  of  the  points  calculated  for  this  event  in  those  figures. 

The  final  data  set,  obtained  several  hours  after  the  first  set  on  August  15, 
presents  evidence  that  either  multiple  stations  or  multiple  interaction  regions  are 
involved  in  energy-dependent  (as  a  function  of  L)  pitch-angle  scattering.  At  the  time 
these  data  were  obtained,  the  OV1-19  satellite  was  located  such  that  electrons  pitch- 
angle  scattered  into  the  drift  loss  cone  west  of  58°  E.  should  have  been  mirroring  above 
the  satellite  and  could  not  have  been  observed.  Particles  precipitated  by  UMS  probably 
could  not  have  been  observed.  However,  particles  were  seen  in  the  drift  loss  cone.  The 
event  is  rather  strange  in  that  the  312  keV  and  250  keV  channels  both  peak  at  L  =  1.66 
and  then  the  250  keV  and  192  keV  channels  both  peak  at  L  =  1.71.  The  two  peaks  in  the 
250-keV  channel  have  different  local  loss-cone  angles  and  trace  back  to  different 
longitudes  of  origin.  The  312  keV  peak  and  the  second  250  keV  peak  trace  back  to  59°  E. 
(which  is  about  the  limit  of  traceback  possible)  while  the  first  250  keV  peak  and  the  192 
keV  peak  trace  back  to  about  72°  E.  Particles  were  also  observed  in  the  139  keV  channel 
but  the  peak  in  L  was  not  observed  because  data  acquisition  was  terminated  prior  to  the 
peak  being  reached.  The  data  are  shown  on  an  East  Longitude  vs.  Universal  Time  graph 
in  Fig.  14.  UMS  was  on  only  during  the  time  period  of  the  origin  of  the  139  keV 
particles.  It  definitely  was  not  on  during  the  times  required  for  precipitation  of  the 
other  energies.  RPS  was  not  on  during  this  time  period. 
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Fig.  14.  Presentation  similar  to  Fig.  1 1,  but  for  data  obtained  on  the  second  pass  of 


August  15,  1969. 


A  SUCCESSFUL  PREDICTION 


Faced  with  the  dilemma  that  electrons  were  present  in  the  drift  loss  cone  with 
the  energy  vs.  L-shell  dependence  characteristic  of  the  VLF  precipitation  model  while  it 
was  certain  that  UMS  was  not  transmitting,  it  was  decided  to  invert  the  model  in  order 
to  obtain  the  frequency  required  to  fit  the  observations. 

The  observation  took  place  one  orbit,  i.e.,  approximately  150  min,  after  the  first 
event  described  above.  The  equatorial  electron  density  for  that  event  is  plotted  as 
squares  in  Fig.  3.  We  assumed  that  the  density  did  not  change  between  the  time  of  the 
two  observations  and  used  the  density  profile  from  the  earlier  pass  with  the  model  to 
compute  the  frequency  of  the  transmitter  required  to  explain  the  energy  vs  L-shell 
dependence  of  the  second  observation.  The  electron  energy,  L-shell,  wave-normal,  and 
particle  pitch  angle  (Fig.  5)  were  fixed.  The  calculations  were  performed  for  the  two 
energies,  190  and  250  keV  which  have  the  best  defined  maxima.  The  results  are  shown  in 
Table  2.  The  frequency  required  is  ~  22.0  to  22.4  kHz. 

Those  frequencies  are  above  the  0  to  20  kHz  range  which  we  had  requested  for  the 
spectrograms  of  the  Byrd  VLF  data.  We  then  requested  a  new  spectrogram  covering  the 
frequency  range  from  20  to  25  kHz.  Four  transmitters  were  operating  continuously  from 
20:34  to  23:15  UT.  Their  identities  and  locations  are  given  in  Table  1. 

The  obvious  candidate  to  explain  the  observations  is  NWC  located  at  North  West 
Cape,  Australia  (114°  E.)  and  operating  at  22.3  kHz. 

As  noted  above,  the  particles  that  were  observed  were  traced  back  to  a  longitude 
between  59  and  72°  E.  There  are  two  possible  explanations  for  this  discrepancy  in 
longitude.  The  first  is  that  the  waves  propagated  westward  in  the  earth-ionosphere  wave 


Table  2 


Parameters  Used  to  Calculate  the  VLF  Transmitter  Frequency 
for  the  Observations  at  22:18  UT  on  August  15,  1969 


Equatorial 

Equatorial 

Wave 

Pitch 

Electron 

Normal 

Energy 

Angle 

L 

Density 

Angle 

Frequency 

190  keV 

28.5° 

1.713 

3850  cm"3 

62.5° 

22.4  kHz 

250 

30.0 

1.662 

4050 

62.5 

22.0 

guide  entering  the  magnetosphere  and  in  fact  interacting  with  the  particles  between  59 
and  72°  E.  Insofar  as  signals  from  the  high  power  transmitters  have  been  observed 
worldwide  by  the  OGO  satellites  [Heyborne,  1966],  this  is  a  possibility.  The  alternative 
explanation  is  depicted  in  Fig.  15.  In  Fig.  15a  the  solid  curve  shows  the  pitch-angle 
distribution  of  stably-trapped  particles.  When  the  interaction  with  a  transmitter  signal 
is  strong  the  precipitated  particles  fill  the  drift  loss  cone  to  the  edge  of  the  bounce  loss 
cone  in  the  appropriate  longitude  range  of  the  transmitter.  The  example  in  Fig.  15a  is 
the  one  typically  observed  where  UMS  fills  the  loss  cone  at  longitudes  west  of  60°  E., 
similar  to  Fig.  12. 

In  Fig.  15b  the  solid  curve  again  represents  the  pitch-angle  distribution  of  the 
stably-trapped  particles.  If  there  is  no  interaction  to  the  west  of  NWC  and  if  the 
interaction  with  the  whistler-mode  waves  is  weak  in  the  sense  that  it  does  not  fill  the 
local  drift  loss  cone  at  NWC  (in  the  example  in  Fig.  15b  it  fills  only  to  24°  which  is 
identical  to  the  case  in  Fig.  15a)  then  the  location  of  the  interaction  region  in  longitude 
cannot  be  determined.  A  priori  there  is  no  compelling  reason  to  believe  that  the  pitch- 
angle  scattering  must  always  fill  the  drift  loss  cone  up  to  the  bounce  loss  cone.  The 
evidence  from  the  resonance  calculations  is  that  NWC  was  the  source  for  the  waves  that 
precipitated  the  particles  observed  at  22:45  UT  on  August  15,  1969.  It  seems  more  likely 
that  weak  pitch-angle  scattering  near  the  longitude  of  the  transmitter  rather  than 
propagation  in  the  waveguide  well  to  the  west  of  the  transmitter  is  the  more  satisfactory 
explanation  for  the  pitch-angle  distributions  observed. 
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Fig.  15.  a)  Idealized  equatorial  pitch-angle  distribution  of  electrons  which  interact 
with  the  atmosphere  at  the  South  Atlantic  Anomaly  (SAA)  and  are  not 
strongly  scattered  by  waves  (solid  line);  modification  to  pitch-angle  distribu¬ 
tion  which  occurs  if  strong  pitch-angle  scattering  occurs  in  the  vicinity  of 
UMS  (dashed  line).  BLC  is  the  bounce  loss  cone  at  the  longitude  of 
observation  and,  for  this  example,  is  assumed  to  be  smaller  than  that  at  the 
longitude  of  UMS  (a  requirement  for  traceback). 

b)  Solid  line  is  the  same  as  in  a),  above.  The  dashed  line  represents  the 
modification  in  pitch-angle  that  would  be  expected  if  relatively  weak  pitch- 
angle  scattering  (i.e.,  insufficient  to  isotropize  the  distribution  at  the 
location  of  the  interaction)  occurred  at  the  longitude  of  the  transmitter  NWC 
(which  has  a  local  bounce  loss  cone  smaller  than  the  bounce  loss  cone  at  the 
longitude  of  observation). 


WAVE-INTENSITY  ESTIMATE 


We  can  estimate  the  pitch-angle  diffusion  coefficient  from  the 
particle  data.  An  estimate  of  the  wave  intensity  required  to  account  for 
the  observed  scattering  can  then  be  obtained. 

The  pitch-angle  diffusion  coefficient  Dao,  is  defined  as 

Daa=  <(6o)2>/2t  (5) 


where  5o  is  the  pitch  angle  change  during  a  single  interaction  and  t  is 
the  average  time  between  interactions  [Ashour-Abdalla,  1972],  In  our 
case  t  =  the  electron  bounce  period.  The  pitch  angle  change  is  averaged 
over  all  initial  phases  between  the  wave  magnetic -field  vector  and  the  com¬ 
ponent  of  the  electron  velocity  vector. 

The  change  in  pitch  angle  dor  that  occurs  during  a  single  interaction 
of  duration  6t  is  [Michael  Schulz,  private  communication,  1980] 


(5 


a)rms  =  y|  bi  ]<iR>  (0  coso 

+  (// C)  [RS  sing  -  COS0]  (C)  coso 


+  [RS  cos9  +  sin8]  (£)  sino|  5t 


+  “  (q/Vm)  bL  \(£/C)  J,  (C) 

v  2 

-  (iR)  '  (C)  cos 0  j  6t 


(6) 
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where  q  is  the  charge  on  the  particle  (for  numerical  evaluation  the  sign 
is  negative  for  electrons),  c  is  the  speed  of  light,  p  is  the  index  of  refrac¬ 
tion  of  the  wave,  p  is  the  momentum  of  the  particle,  R  =  Ex/E^  is  the 
polarization  of  the  wave- field,  S  =  Ez/Ex,  9  is  the  wave-normal  angle  and 
a  is  the  particle  pitch  angle  (both  9  and  a  are  measured  from  the  positive 
direction  of  the  geomagnetic  field;  thus  cos9  in  our  case  is  negative), 
is  the  magnetic  field  intensity  of  a  monochromatic  wave,  and  /  is  the  order 
of  the  resonance.  I  -  +1  satisfies  the  particle  observation.  C  =  k^p^/mil 
where  is  the  component  of  the  wave  vector  perpendicular  to  the  geo¬ 
magnetic  field  B,  p^  is  the  component  of  the  particle  momentum  perpen¬ 
dicular  to  B,  m  is  the  relativistic  mass  of  the  particle  and  ft  is  the  gyro- 
frequency  of  the  particle.  We  use  the  sign  conventions  of  Ratcliffe  [1959] 
where  the  sign  offtis  negative  for  electrons.  For  the  whistler  mode 

iR  «  sgn  [(fto/<j>)  cos9] 

2 

i  (aj  /ou )  ft  sin  9 

g  % - E - 1_2 _ 

[1  -  (ajp/au)2]  (|ftQ  cos  9|  -  cu) 

where  ftQ  is  the  gyrofrequency  for  a  nonrelativistic  electron,  uu/Zn  is  the 
wave  frequency  and  is  the  plasma  frequency. 

Eq.  6  is  a  generalization  for  nonzero  wave-normal  angle  for  Eq. 

2.46  in  Schulz  and  Lanzerotti  [1974].  The  computation  assumes  a  linear 
interaction.  Inan  et  al,  [1978]  has  computed  the  non-linear  pitch-angle 
scattering  of  energetic  electrons  by  coherent  VLF  waves  propagating  with 
zero  wave-normal  angle. 
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The  wave  field  required  to  scatter  a  particle  can  be  calculated  from 


Eqs.  5  and  6  with  appropriate  estimates  for  and  fit. 

The  diffusion  time  required  for  the  pitch  angle  to  scatter  through 
a  total  change  of  Aat  is 

td  =  <(Aa)2>/2  Daa  (7) 

The  observed  change  in  Aa  is  about  3  deg  or  0.05  radians.  A  250 
keV  electron  with  a  pitch  angle  of  28  deg  has  a  bounce  period  of  0.  2  s. 

It  drifts  across  the  30  deg  longitude  sector  of  the  interaction  region  in 
about  300  s.  Taking  =  300  s  and  Aa  =  0.  05  we  obtain  from  Eq.  7 
Daa  =  1.39  y  10"6  s'1. 

fit  is  related  to  the  intrinsic  bandwidth  owing  to  the  inhomogeneity 
in  the  magnetic  field  f Schulz,  1972].  For  equatorial  cyclotron  resonance 
the  intrinsic  bandwidth  isAcu/2Tr,  where  is  given  [M.  Schulz,  1974,  1975] 
by 


| it  <£/2|^  [1  -  (v  /k  •  v  )  cos  9] 


-§ 


(8) 


where 


2  3  2  2 

(3c/uu  La)  (n/v  -  w)  ['jo  (v  +  1  -  sec  a)  +  Q(2/Y  -  vcosG  -  sec  orcosB)] 


and 


[1  -  (v  /k  •  v  )  cos  9]  =  1  +  (n/ Y  -  m)  ficos  9/2tt  ( ft  cos  9  -  c u ) 
ii  g 

for  the  case  of  electron-cyclotron  resonance  for  a  relativistic  electron 
with  a  whistler-mode  wave  with  a  wave-normal  angle  of  9  with  respect 
to  B.  The  cold-plasma  density  is  taken  as  proportional  to  Bv.  The 
"gyrofrequency  model"  used  in  this  paper  corresponds  to  v  =  1.  The  factor 
in  brackets  on  the  right  hand  side  of  Eq.  8  is  a  correction  to  Eq.  4  of 
Schulz  [1972]. 

To  a  stationary  observer  the  length  of  a  wave  train  that  can  interact 
with  the  particle  is  St  [1  -  (v|(/k  •  v  )  cos  9]  in  time.  The  interaction  time 
St  is  related  to  the  bandwidth  A<ju/2tt  by 

Au)/2tt=  (St)  1  [1  -  ( V|  /k  •  Vg)  cos  9  ]~  1  (9) 

For  the  parameter  values  given  in  Table  3  the  interaction  time  5t 
obtained  from  Eq.  9  using  Eq.  8  to  compute  Am  is  2.4  ms.  The  wave 
magnetic- field  intensity  computed  using  Eqs.  5  and  6  is  then  3  x  10"  ^  T 
or  3  mV. 

We  can  estimate  the  magnetic  field  intensity  of  the  waves  at  the 
equator  using  the  simple  analytic  method  outlined  by  Heyborne  [1966]. 

For  an  input  magnetic  latitude  of  30°  the  nighttime  D- region  loss  is  5  db, 
-and  the  divergence  loss  is  7  dB.  Assuming  a  radiated  power  for  UMS  of 
300  kW  [ Watt,  1967],  one  obtains  at  L  =  1.8  a  field  intensity  of  4  mV  for 
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Table  3 


Parameters  used  in  the  pitch-angle 
diffusion  calculation 


Equatorial  Electron  Density,  ng 

3756  cm 

L-shell  parameter,  L 

= 

1.  734 

Electron  Pitch  Angle,  a 

= 

28  deg 

Wave  Frequency,  a>/2ir 

= 

17.  1  kHz 

Wave  Normal  Angle,  0 

= 

1 52. 4  deg 

Energy  of  Resonant  Electrons,  T 

= 

250  keV 

Geomagnetic  Field  Intensity,  B 

s 

5.  75  X  10-' 

Electron  Gyrofrequency,  £2q 

= 

1. 01  X  106 

v/c,  p 

= 

0. 7409 

Electron  Bounce  Period, 

= 

0.  2  s 

Diffusion  Coefficient,  D 

aa 

= 

1. 32  x  10" 

Intrinsic  Bandwidth,  Aco/2tt 

= 

93  Hz 

Wave  Index  of  Refraction,  fj, 

= 

17.  6 

Wave  Group  Velocity,  v 

© 

= 

2.  63  x  107 
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an  all  land  path.  However,  in  the  daytime,  the  D- region  attenuation  is 
35  dB,  which  reduces  the  wave  intensity  to  0.  1  milligamma. 


CONCLUSIONS 


Electrons  observed  in  the  drift  loss  cone  by  the  OV1-19  satellite 
between  L  =  1.  5  and  1. 8  were  pitch-angle  scattered  by  whistler -mode 
waves  from  high-powered  ground-based  VLF  transmitters. 

During  three  time  periods  on  August  14  and  15,  1969  it  has  been 
shown  that  radiation  from  the  VLF  transmitters  UMS  and  NW C  interacting 
with  electrons  and  thereby  precipitating  them  into  the  drift  and  bounce 
cones,  explains  in  a  consistent  manner  certain  particle  observations. 

In  two  cases  the  interaction  with  UMS  was  sufficiently  strong  to  fill  the 
drift  loss  cone  to  the  edge  of  the  local  bounce  loss  cone  at  the  longitude 
of  the  interaction.  In  the  other  case  the  interaction  with  radiation  from 
NWC  did  not  completely  fill  the  drift  loss  cone. 
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LABORATORY  OPERATIONS 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  pKy sics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
Aside;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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